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Abstract 

In these lectures we discuss how the Painleve equations can be written in terms 
of entire functions, and then in the Hirota bilinear (or multilinear) form. Hirota's 
method, which has been so useful in soliton theory, is reviewed and connections from 
soliton equations to Painleve equations through similarity reductions are discussed 
from this point of view. In the main part we discuss how singularity structure of 
■ the solutions and formal integration of the Painleve equations can be used to find 

a representation in terms of entire functions. Sometimes the final result is a pair of 
Hirota bilinear equations, but for Pyi we need also a quadrilinear expression. The 
use of discrete versions of Painleve equations is also discussed briefly. It turns out 
that with discrete equations one gets better information on the singularities, which 
^ \ can then be represented in terms of functions with a simple zero. 
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1 Introduction 



Hirota's bilinear method has turned out to be very efficient in constructing multisoliton 
solutions to integrable evolution equations. But since Painleve equations do not have 
soliton solutions, why should we care about writing them in the Hirota bilinear form? In 
these lectures we will show that this method is relevant even for integrable ODE's. 
The fundamental idea behind Hirota's direct method is the following: 

Change into new variables in which the solutions have the simplest form. 

This transformation will change at least the dependent variable, and may sometimes 
be rather complicated. For solitons the nicest possible form is the one where the soliton 
solution is given as a polynomial of exponentials with exponents linear in the independent 
variables (see. Sec. 2). 

Painleve equations do not usually have solutions that can be written as polynomials 
of exponentials, and although there are other special solutions (rational or solutions made 
of special functions) to which Hirota's method is relevant, there are other more general 
reason that lead to the same. Indeed, the idea of solutions being "as nice as possible" 
can be extended to ODE's: We can demand that the solutions are expressed in terms of 
entire functions. This is not a new idea, it was studied by Painleve himself in |l| |2[]. In 
his Acta Mathematica paper of 1902 he writes [[|, p. 14: 

Puisque les integrates y(x) des equations precedentes [Pj, P/j, and Pui in 
the present notation] sont des fonctions meromorphes dans tout le plan, il 
est bien evident qu'elles sont representables par le quotient de deux fonc- 
tions entieres; mais ce qu'il importe de remarquer c'est qu'on peut choisir ces 
fonctions entieres de maniere qu'elles verifient une equation different ielle tres 
simple du 3 e ordre. 

It should not be surprising that Painleve 's explicit results for Pi, Pu, and Pni in are 
in the bilinear form. More recently solutions to the Painleve equations in terms of entire 
functions were considered by Lukashevich M, and in this school K. Okamoto will give still 
another method of bilinearizing the Painleve equations. 

The outline of these lectures is the following. In Sec. 2 we will introduce Hirota's direct 
(bilinear) method by discussing the soliton solutions of the Korteweg-de Vries equation. 
In this case "niceness" is obvious, because the explicit soliton solutions have the simplest 
possible form. In the subsequent sections we will write the Painleve equations in terms of 
entire functions, using three methods. First in Sec. 3 we use the fact that many Painleve 
equations can be obtained by similarity reductions from soliton equations with already 
known bilinear forms. In Sec. 4, which is the main part, we discuss how quadratic and 
quartic forms can be derived by studying the singularity structure of the solution and 
then writing the equations in terms of entire functions. (For a previous study along these 
lines, see ||.) In Sec. 5 we discuss briefly how discrete Painleve equations (c.f. the talk 
of A. Ramani) can be used as starting points, because somehow the discrete formulation 
is more sensitive to the singularity structure. 
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Finally in this introduction, let us list the equations under discussion: 

Pj : y" = Qy 2 + z, (1) 
P n : y" = 2y 3 + xy + a, (2) 

Pin : y" = - y' 2 - - y' + y 3 + - (ay 2 +P)-K (3) 

y z z y 

eP m : u »=± u '2 + e ^ u 3 + e z(au 2 + f3)-e 2x -, (4) 
u u 

Piv : y" = ^ y' 2 + ^y 3 + W + 2{z 2 -a)y + (3 -, (5) 
2y 2 y 

Pv. y"=(^ + ^-]y' 2 --y' 
\2y y-lj z 



y(y - l) 2 Ay -I) 2 y c y(y + 1) 
z z z z y z y — 1 



ePy : u" = \{- + - it' 2 

2 \M M — 1/ 

it „u — 1 



- ( a^— + p- + -?e x u(u - 1) + 5e 2x u(u - l)(2u - l)\ , (7) 

V u — 1 u J 

2 Vj/ y-i U z-l y-zj 



y{y-l){y-z) ( B ± z-l Mz-1) 

z \z-\y \ "^V 7 (v-i) 2 {y _ z)2 



eP VI : u»=U 1 - + ^ + ^-)u*-e*(^ + 1 



2 \u u — 1 u — e x J \e x — 1 m — e x 



u(u — l)(u — e x ) 



e x /3 (e* — 1)7 e^e* — 1)7 
w 2 (w — l) 2 (w — e x ) 2 



(9) 

The exponential versions are obtained by y(z) = u(x) for Pjjj and Pyi, and y(^) = 
for Py, where z = e x , and the primes of u stand for differentiation with respect to x. 

2 Hirota's bilinear method for soliton equations 

Here we will briefly discuss Hirota's method || for constructing multisoliton solutions to 
integrable equations (for a review, see e.g. 0, 

2.1 Definitions 

The first step in the construction is to transform the equation into the Hirota form. As 
an example let us consider the Korteweg - de Vries (KdV) equation 

Uxxx + QuUx + Ut = 0. (10) 

Let us introduce the dependent variable transformation 

u = 2d 2 JogF, (11) 
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(we will see below that the new function F is regular and simple for soliton solutions) 
and then one can write (|1(]) in the following quadratic form (after one integration): 

F XXXX F — 4F XXX F X + 3F% X + F xt F — F x F t = 0. (12) 

This does not look simpler than ( [To|) but one can write it in a condensed form using the 
Hirota D operator: 

(D 4 X + D x D t )F ■ F = 0, (13) 
where D is a kind of antisymmetric derivative, 

D n J ■ 9 = (d Xl - d X2 ) n f{xiHx2)\x 2=xl= x- (14) 
The minus sign, which differentiates D from Leibnitz' rule, is crucial. We have 

T~\2 II 12 t~\4 llll a I III i o H2 

L) X U ■ U = UU — U , D x = UU — Auu + 6U , 

etc. In the context of ODE's it is worth recalling that Borel and Chazy arrived to these 
expressions by invariance theory 0, [UJ, and observed that they yield equations whose 
solutions are entire. For later soliton computations note that P(D)e px ■ e qx = P(p — 
q)e {p+q)x . 



2.2 Multisoliton solutions 

The KdV-equation is the proto-typical representative of the class 

P(D x ,D y ,...)F-F = 0, P(0) = 0, (15) 

for which the multisoliton solutions are indeed simple in terms of F, as opposed to u. 
The general method of construction multisoliton-solutions is by considering the formal 
expansion 

F = l + e / 1 + e 2 /2 + e 3 /3 + ---, (16) 

where e is the expansion parameter, and truncating this at some order. The vacuum or 
zero-soliton solution (OSS) is given by F — 1. For the one-soliton solution (1SS) only 
one term is needed: it is easy to see, that due to the antisymmetry in ( |I4"D F\ = 1 + e n 
(rj = p ■ x) is a solution of (|T5|) , if the parameters p satisfy a dispersion relation P(p) = 0. 
Fi is the one-soliton solution (1SS) and substitution to ([□]) yields the standard result for 
u. 

The two-soliton solution (2SS) for (|l5l) is obtained from the truncation F 2 = I + /1 + /2, 
where f\ = e m + e m . In order to fix f 2 we note that when we stay in the comoving frame 
of one soliton while the other one goes to ±00 (that is when the other 77 approaches ±00) 
we should get the 1SS again. This means that we should try 

F 2 = l + e^ +e r » +A 12 e^\ (17) 

and substituting this into fll5|) and using the dispersion relation for the parameters Pi we 
find that the equation is satisfied, if the "phase factor" is given by 
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An important point to observe is that the above works for any polynomial P. In fact 
there are still other classes of equations for which the generic form has 2SS, but it should 
be noted that the existence of 2SS does not imply integrability. 

Although 2SS can be constructed for the whole class, 3SS work only for certain equa- 
tions, namely for the integrable ones. It turns out that the ansatz for a possible 3SS is 
fixed by the 2SS: the only one compatible with the 2SS ( |17|) is 

F = 1 + e m + e m + e 773 + A 12 e m+m + A 13 e M * + A 23 e m+ri3 

+A 12 A 13 A 23 e^ + ^, (19) 

where rji = Pi ■ x + 77°, and the parameters pi satisfy the dispersion relation P(pi) = and 
Aij are given by fl!8| ). This form is dictated by the requirement that as one of the solitons 
goes to infinity (i.e., the corresponding rj approaches ±00) the other two should form a 
2SS (0). 



When the ansatz (|L9|) is substituted into Q15D one obtains the condition 
PfaiPi + &2P2 + <? 3 p 3 ) 

<7,=±1 

xP(cr 1 p 1 - a 2 p 2 )P(a 2 p 2 - a 3 p 3 )P(a 3 p 3 - cripi) = 0, (20) 

on the manifold defined by the dispersion relations P(pi) = 0. Since the ansatz was 
completely fixed there are no free coefficients, and since our principle is that we should 
be able to combine any three solitons into a 3SS, we cannot impose any new conditions 
on the parameters Pi either. Thus the condition is on the equation. 

The three-soliton condition (|20| ) can be used to search for integrable equations within 
the class (p~5|) |TT| . Note that in this kind of search there are no initial assumptions about 
the number of independent variables and no preferred time. [This is in contrast with 
searches assuming a structure like u t = F(u,u x ,u xx , ■■■)■} The (nontrivial) results of this 



search follows UTTI: 



(D x -AD x D t + 3D 2 y )F-F = 0, (21) 

(D x D t + aDl + D t D y )F.F = 0, (22) 

(D x -D x Df + aD 2 x + bD x D t + cD 2 )F-F = 0, (23) 

(D x + 5D' x D t -5D 2 + D x Dy)F-F = 0. (24) 

Three of these were known before, Q2~H) is the Kadomtsev-Petviashvili (KP) equation, 
(|22| ) is the Hirota-Satsuma-Ito equation, and (p4f) the Sawada-Kotera-Ramani equation. 
Equation (|2"3"| ) is the only new equation and it is obvious that this equation could not have 
been found by any ansatz assuming simple t dependence. All of these equations have also 
4SS and pass the Painleve test | ]12| . 

Similar analysis of 2SS's and 3SS's have been performed on other types of bilinear 



equations (mKdV|13j, sG|J, nlSpf, and BOg 



2.3 Gauge invariance and generalization to multilinearity 

We have so far discussed Hirota's method only from the soliton point of view, but it has 
been found useful in other approaches as well. In particular the r-functions (=F above) 



have been essential in the Kyoto school approach to integrable PDE's [[16] . 
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One recent important observation is that Hirota forms are intimately related to gauge 
invariance. It is easy to show that if F,G, . . . solve some bilinear equations, so do 
e ax F,e ax G, . . .. But the reverse is true as well [17]: If some quadratic expression is 
gauge invariant, then all derivatives must appear as Hirota derivatives. The proof is sim- 
ple. Consider the quadratic homogeneous combination A n (f, g) := Yh=o c * (®x /) (®x~ l flO ■ 
From the gauge invariance A n (e ax f , e ax g) = e 2ax A n (f,g) we can solve for the constants 
c, and find that q = (— l) l (^co, so that A n can indeed be written in terms of bilinear 
derivatives D: A n (f, g) = c D%f ■ g. 

This gauge principle can be applied to higher multilinear expressions |17]]. For the 
cubic case one finds that for gauge invariant expressions the derivatives appear through 

T = d 1 + jd 2 + fd 3 , T* = d x + fd 2 + jd 3 , (25) 

where the subscript indicates on which factor the derivative operates, and j = e 2 " r/ ' 3 . The 
multilinear generalization is 

n-1 

M™ = J2 e 27Tikm/n d k+1 , where < m < n. 

k=0 

One can now search for integrable equations from the class 

P(T, T*)F ■ F ■ F = 0, (26) 
and new equations have been found in ||18|| , for example a generalization of the KP equation 
(T*T* + 8 T^T y T* + 27 Tjj - 36 T*T t )F ■ F ■ F = 0, (27) 
or in the non-linearized form obtained with F = e 9 , 

gxxxxy ~\~ 8g X xygxx ^gxygxxx ^9yyy ^9xxt 0. (28) 



3 Bilinear forms and similarity reduction 

Similarity reductions of PDE's to ODE's of Painleve type are very important theoreti- 



cally, in fact the ARS conjecture |nj] states that if an integrable PDE is reduced to an 
ODE, the ODE should be of Painleve type. We will now follow this reduction path, 
but with a different purpose: Since the bilinear formalism has been so useful and is well 
known for soliton equations, we will use them as starting points and then apply simi- 
larity reductions in order to derive bilinear forms for ODE's. We will not consider here 
all possible similarity reductions to Painleve equations, but just some typical cases with 
direct reduction. [In many case the connection to Painleve equations goes through rather 
complicated (differential) transformations, which probably does not help in the present 
objective of getting bilinear forms.] For further references about similarity reductions, see 
|3, Sec. 6.5.15 and 7.2. 



For bilinear variables the similarity reduction is always assumed to be of the form 
F(x,t) = (p(z)e a ( x,t \ ... where the exponents are to be determined so that the bilinear 
equation is in terms of z only. 
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3.1 Pj 

To get a similarity reduction to Pj let us consider the KdV equation (|i~0l). If one substitutes 
into it the ansatz pi 

u(x,t)=2t-2y(z), z = x-6t 2 , (29) 
then one gets for y(z) the equation 

y"' = 12yy' + l, (30) 

which integrates to Pj ([[]). 

As was shown before, the bilinearization of KdV proceeds through the dependent 
variable transformation (11) so that we now have the relation 

t-y(x-6t 2 ) =d 2 x \ogF. (31) 

This suggest that for Pj we should introduce a new dependent variable cf> by 

</ = -(log </>)", (32) 

and from ( |3"1~| , |3"2"| ) we find that the similarity reduction for the bilinear dependent variable 
corresponding to fl29"D should be 

F = (f)(x — Q t 2 )e^ tx2+a{t)x+m . (33) 

(Note the free functions a and b, on which we have no information at the moment.) When 
this is substituted into ( Jl3| ) we obtain something that is a function of z alone, if we choose 
a(t) = —At 3 (b drops out). The result is then 

{D\ + 2z)4> -0 = 0, (34) 

which is the standard bilinear form for Pj. The notable feature in the above process is 
the necessity of the gauge factor, in this case e^ tx ^ x ~ 8t \ 
We could also start from the Boussinesq equation 

u xxxx "I - 3(it ) xx + u xx Utt — 0, (35) 

and using similarity reduction u = —2y(x — t) |fS| we immediately obtain y"" = Q(y 2 )" 
which can be integrated twice to yield ([!]) with suitable integration constants. Equation 
fl35| ) can be bilinearized as KdV with (|TTD, which yields (after two x integrations) 

(Dt + D 2 X - D 2 )F ■ F = 0. (36) 

The similarity reduction for F should now be of the form 

F = <j>(x - t)e xa{t)+Kt \ (37) 

and indeed the bilinear form (|3^) follows, with z = x — t, if we use the gauge e~^ t2x+ ^ tS . 
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3.2 P, 



ii 



The second Painleve equation can be obtained by a similarity reduction from the mKdV 
equation 

u xxx - 6u 2 u x + u t = 0, (38) 



by the reduction ansatz [22 



u = (3t)- 1/3 y(z), z = x/(3t) 1/3 . (39) 
This yields for y the equation 

y"' = 6y 2 y' + zy' + y, (40) 

which can be integrated to @. 

There are two ways to bilinearize mKdV. In the conventional approach we have to use 
the potential form, as for KdV. Thus let us introduce v by u = d x v, substitute this into 
and integrate the result with respect to x, this yields 



v xxx - 2(v x f + v t = 0. (41) 

(The integration constant can be absorbed into v, since it is defined up to an additional 
function of t.) The bilinearizing dependent variable transformation is 

v = log j (42) 

and substitution into fl4lD yields 

- FG[(D 3 X + D t )F ■ G] + 3[(D 2 X )F ■ G][D X F ■ G\ = 0. (43) 

At this point we have one equation for two functions, so in principle we can introduce extra 
conditions for them. Recall that F and G are defined only up to a common multiplicative 
factor, so this is the origin of the freedom we now have. For soliton solutions it turns out 
that the best way to fix this factor is to demand F)\F ■ G — 0, then we get the bilinear 
form 

/ (D x + D t )F ■ G = 0, 

\ D 2 X F-G = 0. l44J 

The 1SS for this class of equations is given by F = 1 + e v , G = 1 — e v with dispersion 
relation given by the odd polynomial. 

Maybe a general comment on equation splitting is in order here. It should be noted 
that for some other kind of solutions the above might not be the best way to split (f§3|). 
The general method is to put D\F ■ G = XFG where A is an arbitrary function, which 
yields the pair 

f (D 3 x + D t -3X)F-G = 0, 
\ {D 2 - X)F ■ G = 0. 

If we now make a gauge change 

F -> e e F, G -> e e G, 



S 



the above equation changes to 

j {Dl + D t - 3(A - 29))F ■ G = 0, 
1 (Dl - (A - 29))F ■ G = 0. 

For a given type of solution (rational, soliton) one needs a specific form of (A — 26), for 
soliton solutions this term should vanish. 

The other bilinearization of (|38|) is obtained by substituting u — g/f directly into it, 
and the result can then be split into an nlS type bilinear equation 



(Dl + D t )f-g = 
Dlf-f + 2g* = 



2, , , o„2 _ n ( 45 ) 



The 1SS of this system is given by / = 1 — e 2v , g = —2pe v . 

The dependent variables of these two forms ([|4],[|5]) are related by 

g = D x G-F, f = GF. (46) 



Let us now see how the above bilinear forms can be used to bilinearize Pn- In the 

A 

dz 



first case with bilinearization through u = d x \og(G/F) the natural ansatz is y — log ^, 
because then 



^ iog ^ =M= (^^ )= (3^^ iog ! = ^ log !' 

(note the partial derivatives) so that we could just try 

G{x, t) = tfj(z), F(x, t) = <f)(z), with z = x/(3t) 1/3 . (47) 
Indeed this works, and we get from ( [44]) 

(48) 



(Dl-zD z )<p-i) = 0, 
Dl<t>-^ = 0. 



In the second case with u — g/f 



9 1 ( , 1 Hz) 



/ (3ty/* yK ' {3t) l / 3 ^z) , 

suggests we should try 

g = a(x, t)^f(z), f = a(x, t) (3*) 1 / 3 $(z), (49) 
and then from (|45| ) we get a bilinear form depending only on z, if we just choose a = 1: 

j (D»-*Z>. + l)*-* = , 

\ J D z 2 $-$ + 2^ 2 = . [M) 

It is easy to check that the substitution \1/ = <jnf), $ = D z (f> ■ tp reduces (|50[) to (|48|) . 
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3.3 Pin 

Special cases of Pm can be obtained by similarity reductions |22|, |23| from the sine-Gordon 
(sG) equation 

u xt = smu. (51) 

The first similarity ansatz is 

u(x,t) — — Hog y(z), z = xt, (52) 
and substitution to (^L|) leads to the the special case 

y"=7» a -ly' + TM-v- < 53 > 

The sG equation (^1|) can be bilinearized using 

u = -2ilog(±^, (54) 
J ~W 



yielding 



, (D x D t - l)g ■ f = , 
\ D x D t (f.f-g.g) = . ^ 

The similarity reductions for /, g should be g(x,t) = <fi(z), f = ip(z) and they yield |23| 

ip + icj) 



V=iPt\\ (56) 



(K 



with 

f (zD 2 z +d z -l)<j ) ^ = , 

\ (^ 2 2 + 9 2 )(0-0-^-^) = o . 

This, however, is not satisfactory, because contains ordinary derivatives. The trick to 
eliminate them is to change the dependent variables by <f>(z) = 0(£)> i'i 2 ) = ^(0> z = e ^ 
because then we get 

r (d|-o0-^ = o , 

\ J D|(0-0-^-^)=O . 
Another similarity ansatz for ( |5"1"D is 

u(x, t) = -2i log w((), C = 2Vxi, (59) 

leading to 

w " = - w ' 2 -- w > + l( w z -l/ w ). (60) 

w z 

This is related to the above as follows: If 0(£), solve fl5*8|) then 

_ ^(21og(C/2))+^(21og(C/2)) 
W ^(21og(C/2))-^(21og(C/2))- 1 j 

Thus we have the same basic bilinear equation ( j58|) corresponding to two different non- 
linear ones. 
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4 Solutions in terms of entire functions 



As was mentioned before, Painleve considered already quite early the question of repre- 
senting the solutions in terms of entire functions |3|. But how could we find such 
entire functions? (Painleve does not give any constructive method, but just the solu- 
tions.) One direct way is by studying the singularities of the solutions and then doing 
some manipulations on them so that their entireness becomes manifest || . 

Here we would like to present an additional aspect to the introduction of the entire 
functions: By choosing these functions properly one can actually integrate the equation 
once. 

Suppose we have an equation of the form 

y" = m/ 2 + /V + 7, (62) 

where a, j3, 7 are functions of z and y, primes stand for derivatives with respect to z. We 
want to integrate it to the form 

I:=Ay' 2 + By> + C-J c Z Dd(, (63) 



i.e. to find functions A, B, C, D of z and y such that 

§ = y"(2Ay' + B)+A y y' 3 + (A z + B y )y" , K „ M , „ v)y , „ M ^ 



dl = y »(2Ay' + B)+ A y y' 3 + (A z + B y )y' 2 + (B z + C y )y> + C z - D 
= {2 Ay' + B){y" - ay' 2 - (3y' - 7) = 0. 



(Here the subscripts stand for partial derivatives.) This immediately yields the set of 
equations 

A y = -2aA, 
A z + By = -2pA-aB, 
B z + C y = -2>yA-pB, 
C z -D = - 7 S. 

In the following, it often turns out that 

f-= e II Ddzdz (66) 

is an entire function, and that the other entire function can be obtained from g := yf. 
With this definition of / we get two equations from the above: 

Q± = Qo S f)"-D{z,g/f) = 0, (67) 



and 



- = (log/)' - f Z Ddz 

Q Jc 

= (log/)' - [A(z, J )(J )' 2 + B(z, J )(f )' + C(z, J)] - d = 0. (68) 

where C\ is a constant. 

Below we will show that for the Painleve equations R defined above is quartic (with 
q a simple quartic polynomial of / and g (no derivatives)) and Qi quadratic in /, g and 
their derivatives. From these two equations further equations can be derived, including 
those in Hirota bilinear form. 
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4.1 Pj 

For Pj the situation is a special and one entire function is enough. Using the above 
method we find that with a = (5 = and 7 = 6y 2 + z one solution to (|65"D is given by 

A = i, 5 = 0, C = -(2y 3 + ^), £> = -</. (69) 

Painleve mentions that / := eff Ddz is entire when / D dz = \y a — 2y 3 — yz + ci, in 
accordance with (|69|). This is easy to prove: Near any singularity the solution y of ([3]) 
behaves as || 

V= , 1 , 2 +O((z-z ) 2 ), (70) 

[Z — Z ) 

so that at that point / [as defined by with D = — y] behaves smoothly: 

/ = (z - Zq) ■ [const + 0(z - z )). (71) 
Then from ( |67|) (using y in place of p//) we get 

y=-(log/)", (72) 
and when this is substituted into Pj we get for / an equation in Hirota's bilinear form 

(Dt + 2z)f-f = 0. (73) 



In this paper we also want keep track of the integration constants. Equation (ff3|) is 
fourth order so it has two additional constants of integration. They are related to the 
gauge invariance under / — ► e a+x/3 f, which is a common property of equations in Hirota 
form (in [17] is was argued that the gauge invariance is the defining property of Hirota 
form.) 

In the following the final results for other Painleve equations cannot be written as one 
quadratic equation but rather as a pair, so let us do it here also. For this purpose we 
take another solution of (p5|) with A, C, D multiplied by 2 of what was given in in (|69|). 
This yields / = exp(— 2 / dz J dz y) = [z — z ) 2 ■ [const + 0(z — z )] which is needed to 
guarantee that g :— yf is also entire. Then from (|67|) 



Q x = f"f - f 2 + 2fg = \D z f ■ f + 2fg = 0, (74) 
and from fl68|) (g = — f 4 ) 

R = (f'g - g'ff - ff - Ag 3 f - 2zgf - Cl / 4 = 0. (75) 
These equations are equivalent to ([I]) in the following sense: 

2{fg-fg')P I = Q 1 + f(R/f)'. (76) 

The pair (f74,75|) is third order, and since two constants of integration are accounted for 
(ci and the overall scale of /, g) only one more constant of integration remains, and in 
this sense this pair represents the once integrated Pj. 
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Further equations can be derived as follows: By considering (R/(f 2 g 2 ))' = and using 
fl?4] ) to simplify the result we get another quadratic equation 

Q2 = 99" ~ 9' 2 + //' + zfg + cj 2 = 0, (77) 

which appears in ||. However, it is not gauge invariant and therefore not expressible in 
Hirota form, furthermore the pair Q\ = 0, Q2 = is not equivalent to ([I]) (one would need 
R = as well). If one instead considers the combination Q 3 := (g 2 Q\ + PQ2 + R)l(fg) 
one obtains 

Q 3 = f'g - 2f'g' + fg" - zf - 2g 2 = D 2 J ■ g - zf - 2g 2 = 0, (78) 

and ([I]) is equivalent to Qi = 0, Q 3 = 0. Furthermore this pair is in the Hirota form, and 
the two integration constants are related to the gauge invariance (/, g) (e a+x @f, e a+x ^g). 

Thus in terms of the entire function / and g Pi can be expressed by one fourth order 
equation in Hirota form (73) or by the third order pair (74,75), or by the pair of second 



order equations in Hirota form ( |7^J78| ). 
4.2 P n 

For Pi 1 the expansion around a singularity is || 

y = ±^—T^(z-z ) + ... (79) 
z — z 

and if we just consider y 2 we get entire functions from 

f;= e -ffy 2 <i^z = z - Zo + ..., g:=yf = ±l + ... (80) 

The integration yields the solution 

A=l, 5 = 0, C = -(y 4 + zy 2 + 2ay), D = -y 2 . (81) 

agreeing with the above. (Painleve considers J D dz = y' 2 — y 4 — zy 2 — 2ay in [l], 
Then from ( |6^ ) we get the equation 

Q x = ff" - f 2 +g 2 = \D 2 J -f + gg = 0, (82) 

and from ( |68"1) (g = — / 4 ) 

R = (f'g - g'f) 2 - ff -g'- zg 2 f - 2agf - c x f = 0. (83) 

(both given by Painleve in ||). Equation ( |76|) holds also for Pn. 

As before, another quadratic equation 0] is obtained from [Rj (fg 2 ))' = 

Q2 = gg" - g' 2 + ff + agf + cj 2 = o. (84) 

(Note that here z is absent.) Instead of this one could consider the gauge invariant (and 
ci independent) expression Q 3 := (g 2 Qi + fQi + R)/(fg), i-e., 

Q 3 = f'g - 2/V + fg" - af - zfg = (D 2 Z - z)f ■ g - af = 0. (85) 
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The Hirota bilinear pair (32,35) is the same as given in ||. The counting of integration 
constants is as before. 

Still another form is obtained if we take / = FG, g = D Z F ■ G corresponding to 
y = F'/F — G'/G 0, which leads to the bilinear form 



{Dl 



D 2 F-G = 0, 
zD ? -a)F-G = 0. 



(86) 



This is fifth order and there are 3 obvious integration constants related to the invariance 
under F — > ae cx F, G — > be cx G. 



4.3 P IU 

For Pin the expansion around a movable singularity is 

y = ± — - + . . . (87) 

Z — Zq ZZo 

and if one consider the combination z(y 2 + -y) = Zq/{z — z ) 2 + 0(1) on finds that 

/: = e -/f/^+f^, g:=y f (88) 

are entire p. 

The term — above suggests that it might be better to work with the exponential 
version (Q) (as was done by Painleve, 0). Then one solution to the integration problem 
is 

A = l/u 2 , B = 0, C = 2e x {(3/u-au) - e 2x (1 / u 2 + u 2 ^ 

D = 2e x (p/u-au) -2e 2x (l/u 2 + u 



2,n u.2 , .2^ (89) 



This corresponds to Painleve's 2£ in p. 15. However this does not directly lead to 
entire functions, and Painleve adds some ad hoc operations, which in fact amount to 
using another solution 



A = l/(4w 2 ), B = — l/(2u), C = \{e x {(3/u-au) - \e 2x {l/u 2 + u 2 )), 

D = —(e x au + e 2x u 2 ). 



(90) 



This leads directly to the desired result: / := and g := uf are entire, and we get 

Q 1 = f'f - f' 2 + ae x fg + e 2x g 2 = \D 2 J ■ f + ae x fg + e 2x g 2 = 0, (91) 
and (q = -4/V) 

R=(fg- g'ff -2fg(fg+ g'f) + fg 2 - 2e x (afg 3 - (3fg) - e 2x (g 4 + f 4 ) - Ac.fg 2 . (92) 

As usual, by considering (R/ (f 2 g 2 ))' = we get another equation, which now happens to 
be in the Hirota form: 

Q 2 = 99" ~ g' 2 - Pe x fg - 5e 2x f 2 = \D 2 x g ■ g - (3e*fg + e 2x f 2 = 0. (93) 
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Thus Pin is equivalent either to the pair (|9l"|j92"l) or (|9T| , |93D , in the first case the system is 



third order and there are two integration constants, the overall scale and Cj, in the second 
case the system is fourth order with two-parameter gauge freedom. 

Note that Pm is invariant under u — > 1/u accompanied with the parameter changes 
a — > —j3, (3 — > —a. This corresponds to / <-> g and we see from the above that it is 
indeed a symmetry of the bilinear equations. Thus it might be said that the zeroes of u 
are as important singularities as its poles, and to handle all of them at the same time one 
could define functions F, G, K, M by || 

G' F' e x M' K' 
f = FG, g = KM, e* M =---, --jf-jf- (M) 

In that case we get four equations for four entire functions, two from the above definitions 

D X G-F = e x KM, D X M ■ K = e x FG, (95) 

and two from Qi, Q 2 

D 2 X F-G = -ae x KM, D 2 X M ■ K = f3e x FG. (96) 

The system is now 6th order with four scale related integration constants: F —>■ abe cx F, G - 
a/be cx G, K =—>■ kme cx K,M =—> k/me cx M. The result (|95|j96|) is more symmetric, but 



involves twice as many dependent variables. Whether it is more useful in practical appli- 
cations depends on the problem. 

4.4 P IV 

For Piv the expansion reads 

y = -z Q + ... (97) 

z- z 

and one finds that 

f:=e -f*J*<V+*y) j g . = yfi (98) 

define entire functions @. 

The integration method again works with the simplest choice. If we use the solution 



A = —, B = 0, C = -\(y 3 + 4zy 2 + 4(z 2 -a)y + 2(3/y), D = -(y 2 + 2zy), (99) 



1 

4? 

agreeing with fl98|) , we get the equations 



Q 1 = ff - f 2 + g 2 + 2zfg = \D 2 J ■ f + g 2 + 2zfg = 0, (100) 



2 2 



and (g = — 4/ 3 

R = (f'g - g'f) 2 - Afgf - g A - Azfg 3 - A{z 2 - a)fg 2 + 2(3 f - 2c, fg = 0, (101) 
and the other quadratic equation is 

Q 2 = g"g - g' 2 + 2gf - 2(3 f + cjg = 0. (102) 



15 



The gauge invariant combination of the above turns out to be trilinear: 

T = T Z T*J -g-g- 2{(3 f + 2{z 2 - a)fg 2 + zg 3 ) = 0. (103) 

(T z T*f ■ g ■ g = gD z g ■ f + \ fD z g ■ g.) One can now show that P IV can be expressed as 
a linear combination of either Q\ and (R/(f 3 g))' or Qi and T, with the usual accounting 
of integration constants. 

At this point we would like to return to the question of gauge trasformations, briefly 
mentioned before. The point is that the function 

j .— g- J dz J dz(y 2 +2zy)+p{z) 

is entire for any fixed polynomial p of z. For example note that Pjy has the polynomial 
solution y = — |z, and then from (log/)" = — y 2 — 2zy we would get / cx e^ z . It would 
clearly be desirable to have polynomial /, g as well. This could be obtained by a proper 
choice of p, see p6|, p. 68 for details. The same problem exists for Pni, see P6fl, p. 90. 



4.5 Pi 



V 



In this case again the nicest results are obtained for a specific form of the equation. 
Computations with the standard form reveal that one should instead consider the equation 
(|7|) obtained from the standard one by y(z) = u(x)/(u(x) — l),z — e x ||. The expansion 
for u is given by 



^ ± i/Vg_ ± ^ + 7 -2fc + 

Z — Zq Aozq 
and using the method of [[J leads one to the entire functions 

/ := e f dx f dx(ye*u+26e^u(u-l) ^ g ;= U f . (105) 

The integration has the corresponding solution 

A = - . 1 rr , B = 0, C = -a—^—-(3- + 1 e x u + 6e 2x u(u-l), nnR , 
2u{u — 1) u — 1 u (106) 

D = 'je x u + 2Se 2x u(u — 1), 



and from (p7|) we get the first equation 

Q, = \D 2 J ■ f - ie x fg - 25e 2x g(f - g) = 0, (107) 

and from © ( e = -2fg(g - /)) 

i? = (fg-fg')i-2fg(g-f)f-c 1 f 2 g(g-f) 

-2 f[ag + (3{g - /)] + 2( 5 - Z)^^/ - 5e 2 *(s - /)] = 0. (108) 

From the derivative (R/(f 2 g 2 ))' = one obtains equation 

Q 2 = «/'<? - - (,/' - 2/3/ 2 + {a + /3- Cl /2)fg = 0, (109) 

and the gauge invariant linear combination of the above is again trilinear, 

T = T x T:(f-lg)-g-g-2afg 2 -2(3f(f-g) 2 - 1 e x g 2 (2f-g)-25e 2x g 2 (f-g) = 0. (110) 

(T x T*(f - \g) ■ g ■ g = gD 2 J ■ g + (//2 - g)D 2 x g ■ g) Furthermore one finds that P v is 
expressible as a linear combination of Qi and (R/(f 2 g(f — g)))' or of Q\ and T. 
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4.6 P VI 

For Pyi the situation is more complicated. In fact the method used in || does not work: 
there are no polynomials of u alone from which entire functions can be built. On the other 
hand Painleve in |24| proposes an expression that is supposed to yield an entire function, 
but this expression involves also u' . Presumably one can search such expressions using 
the expansion around the singularity, but we will here take a different route. 

It turns out that Lukashevich [|J obtained some quadratic and quartic expressions for 
Pyi as well, but we could not verify the precise forms given in ||]. Using these results as 
a guide we searched for two quadratic and a quartic expression with similar properties as 
before, using This resulted in 

Qi := (e« - If [f'f - f' 2 ) + (e* - l)fg' + 2ag(g -/)-(<* + c x ){e* - l)fg, (111) 

Q 2 : = e-(e^-l) 2 (^-^ 2 ) + (e^-l)/^ + /3eV(^-/)-(/3-5-7 + ci)(e ;c -l)^, (112) 
and 

R := (e x -l) 2 (f'g-fg') 2 -2(e x -l)fg(f-g)(e x f'-g') 

-2ag 2 (f - g)(e x f - g) + 2^f\f - g){e* f - g) - 2 7 (e* - l)fg(e x f - g) 



-2Se x (e x - l)fg{f - g) + 2 Cl (e x - l)fg(f - g)(e x f - g). 

The relationships between these expressions and the Pyi equation are as follows: 

2(e*-l) 2 fg(f-g)(e*f-g)P VI 

= 2{-g 2 Q x + e*f 2 Q 2 )(f - g)(e x f - g) + (e*f 2 - g 2 )R, 

and 

MP VI = 2e x (Q 1 - Q 2 ) + (e x f 2 - g 2 ){e x - l) 2 (R/U)', 
where 

U := ( e *-l)fg(f-g)(e x f-g), 

M : = -2(e x -l) 3 f 2 [(e*-l)(f'g-fg')(e x f 2 -g 2 )-e*fg(f-g) 2 ]/U 

(Note a spurious solution: M vanishes if u solves (e x — l)u'(u 2 — e x ) = e x u(u — l) 2 .) 
Finally we have a relation between Qx, Q 2 and R as 



(113) 



(114) 



(115) 



where 



B 
C 
V 



BQ 1 + CQ 2 + (e x - l)fgV(R/V)' = 0, 



-2g 2 i(e x -l)(f'g-fg')-e x f(f-g)], 
2f 2 [ e -(e x -l)(f'g-fg')-e x f(f-g)) 

\2 r2„2 



;ii6) 



: (e^-l) 2 /V 

As far as gauge invariant expressions are concerned, one finds that 

XQi + YQ 2 + ZR 
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is gauge independent whenever 

X + Y = 2Z(f-g)(e*f-g), 

and then the c\ terms vanish as well. One possibility is to take 

Q x -Q 2 = l(e B -l)\Dlf.f-e-*Dlg-g) + (e B -l)D x g-f 

+2(ag - (3e x f)(g -/)-( a -/? + 7 + 6)fg(e x - 1) = 0, (117) 

which is bilinear. For the other expression we could not get any simplification so it is 
quadrilinear, for example: 

{e x -lf{f-9){e x f-g)Dlf-f 

+ (e x - 1) 2 (DJ ■ g) 2 - 2e x (e x - - g)(D z f ■ g) 

-2ag(f - g)(e x f - g)((e* - l)f - g) + 2f3e x >f(f - g)(e x f - g) 

-2j(e x - l)fg{e*f - g) - 25e x (e x - l)fg(f - g) = 0. (118) 

Let us now return to the integration procedure and try to understand why the straight- 
forward procedure failed. Substituting g = uf into (|113[) shows that 

U (log/)'- {Au' 2 +Bu' + C]- Cl , (119) 



-2{e*-l)fg(f-g){e?f-g) 
where 



e x — 1 „ —1 ^ mi e x 3 7 e x 5 

A = — r -T7 Z\i B = -, C = — - — - + — - — t^ + —- + 



2u(u — l)(u — e x ) ' u — e x ^ e x — 1 (e x — l)u u — 1 u — e x 

(120) 

This expression is in fact in p4}| , (eq. (3), m = Au' 2 + Bu' + C from above), and Painleve 
states that e^ m has no singularities, apart from the fixed ones (for z they are at 0,1 and 
00, for x{= log z) at —00, 0, 00). 

The integration procedure therefore works as before up to this point, and the problem 
is in D of (|S3"D, it is no longer a function of u only. Indeed, if one writes again 

I := Au' 2 + Bu + C - A (121) 

and identifies A with (log/)' then from Qi we get 

(e x - 1)A' + u {A-a- c x ) + 2a ^~ 1} = 0. (122) 



(Similar expressions were considered in 

For P VI the situation has then turned out to be quite different from the others, as 
might have been expected. Nevertheless, even in this case the final result can be written 
in multilinear form, in this case we need one bilinear ( 117|) and one quadrilinear ( |118| ) 
expression. 
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5 Discrete Painleve 



At the moment the most interesting developments in the field of integrable systems seem 
to take place in the area of integrable difference equations. Most properties of continuous 
integrable systems can be extended to the discrete case, e.g., Lax pairs, existence of 
solitons (for partial difference equations) and the Painleve test. [For an overview see the 
lectures of Nijhoff and Ramani.] 

In order to define discrete Painleve equations one should have a definition of discrete 



Painleve property. Grammaticos, Ramani and Papageorgiou f27[] have proposed that sin- 
gularity confinement is the proper discrete analogue of the Painleve property. Singularity 
confinement means that if a mapping leads to singularity, then after a finite number of 
steps one should get again out of it and this should take place without essential loss of 
information. Singularity confinement has subsequently been used to generate discrete 
forms of Painleve equations, in fact several families of them. (Before calling some differ- 
ence equation a discrete version of a differential equation, one must verify at least that 
its continuum limit is the original continuous equation, but the continuous and discrete 
equations should share some other properties as well.) 

The bilinear approach has a natural analogue in the discrete case. It is best stated 
using the gauge principle: If the expression is homogeneous in the dependent variables 
F,G, . . . and invariant under F(n) — > F(n)e pn , G(n) — > G(n)e pn , . . . then we say it is in 
Hirota form. 



As an example let us consider d-P/[30j. One version is 



W + w + w = ha, (123) 

w 

where w = w(n+ 1), w = w(n),w = w(n — l) and z = an + (3. If in this equation one hits 
a singularity, it is by first arriving somehow to w{k) = 0, and a closer study indicates that 
the sequence of special w values are {0, oo, oo, 0}, after which regular values are again 
obtained. This pattern of w values is obtained from 

, , F(n + 2)F(n-l) , . 

win) = ;T ; r - ; , 124 



if F has a simple zero F(k — 1) = 0. The expression ( |1 24j ) is homogeneous and gauge 



invariant, and if one substitutes it to the discrete derivative of (11231) one arrives to 



F(n + 3)F(n - l)F(n - 2) - F(n + 2)F(n + l)F(n - 3) = 

z(n)F(n + l) 2 F{n - 2) - z(n - l)F(n + 2)F(n - l) 2 . (125) 

This is the trilinear version of d-P/. 

If one now applies the continuous limit to the above, the previous results are obtained 
33: If a = 6, w = 1 + e 2 y, z = -3 + e 4 (, C = ne, one finds that y = 2(logF) c? and flT25|) 
becomes the z derivative of fl73[) divided by f 2 . 

The most important aspect of the above is the way the complicated singularity pattern 
of w is obtained from a simple zero of F. This is the discrete analogue of expressing the 
original solution in terms of entire functions. In the discrete case the process is much 
more clear and this is an indication that discrete systems are more fundamental. 
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If the system has several singularity patterns we need more functions to handle them, 
in general the number of singularity patterns is the same as the number of entire functions. 
In this idea was followed to its logical conclusion: for d-Pyi the authors obtained a 
set of bilinear equations involving 8 functions. We will not repeat all of their results here, 
just some illustrative examples. 

One version of d-Pu is 

- , zw + a 

w+w — - -. (126) 

1 — ur 

It has two singularity patterns, {— l,oo, 1} and {l,oo,— 1}. The entrance to the first 
pattern and exit from the second is described by 

/ s F(n + l)G(n - 1) , , 

»<"> = -»+ hwk) (127) 

while for the remaining part we would get 

Equating these two expressions we get the first equation 

F(n + l)G(n - 1) + F(n - l)G(n + 1) - 2F(n)G(n) = 0, (129) 



whose continuous limit yields the first equation of (|86|) . The second equation is obtained 
from (|126|) . In this case we have two singularity patterns and functions, and the structure 
of the patterns determines one equation between the functions. 

For the higher discrete Painleve equations the singularity patterns sometimes deter- 
mine everything. More precisely, the singularity patterns suggest different ways of writing 
the dependent variable in terms of functions with simple zeroes, and comparing these ex- 
pressions one gets enough equations. As an example let us consider d-Pm 

cd(w — az)(w — bz) 

ww = -) y } 130 

[w — c)(w — a) 

where z = A™ (corresponding to change of variables z = e x in the continuous form) the 
singularity patterns are {c, oo, d}, {d, oo, c}, {az, 0, bz}, and {bz, 0, az}. This suggest the 
representations |30| 

F(n+l)G(n-l)\ J f i F(n-l)G(n + l)\ HK 

(131) 
(132) 




F(n)G(n) J \ F{n)G{n) J FG' 

H(n + l)K(n-l)\ 1 / H(n - l)K(n+ 1)\ FG 



H(n)K(n) J bz \ H(n)K(n) J HK' 

Equating these expressions and taking suitable sums and differences yields four bilinear 
equations in a nice symmetric form 

F(n + l)G(n-l)+F(n-l)G(n + l)-2F(n)G(n) = - (- + H(n)K(n), 
F(n + l)G{n-l)-F(n-l)G{n + l) = ( - - -) H(n)K(n), 



,d c 

H(n + l)K(n-l) + H(n-l)K(n+l)-2H(n)K(n) = -z(a + b)F(n)G(n) 
H(n + l)K(n - 1) - H(n - l)K{n + 1) = z(b - a)F{n)G{n). 
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The continuous limit is obtained with z = e — e , a = e — a e , 6 = — e — 6 e , c — 
1/e + Co, d = — 1/e + d and yields 

D X F-G = -HK, 

D X H ■ K = -e 2x FG, 

D 2 X F-G = (c + d )HK, 

D 2 X H-K = (a + b )e 2x FG. 

One can now verify, that if one uses substitution u = — e~ x '^ log(F/G) in (|) the result 
vanishes due to the above equations. 

The important point in the above construction is that the singularity patterns deter- 
mined all of the bilinear equations. The Painleve equation is then just a way to represent 
the singularity patterns through one function and its equation. This situation continues 
with some of the higher discrete Painleve equations. 

The approach of taking the discrete singularity patterns seriously and using them to 



derive bilinear forms |30| is systematic and powerful. In this way the Hirota bilinearization 
of Pyi was first obtained. The drawback is in the proliferation of dependent variables: 
for Pyi 8 functions are needed, one for each singularity type (the singularities are at 
y = 0,oo,l, and z, with the next term in the expansion having ± sign). In the bilin- 
ear+quadrilinear form ( |117| , |118| ) we manage with two functions, and this is enough also 
for Okamoto's method (truly bilinear but using also ordinary derivatives and therefore 
not gauge invariant). Which form is best will then depend on the practical problem on 
hand, and it is useful to keep all alternatives in mind. 
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